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Abstract 

It is tempting to raise the issue of (metric) chaos in general relativity since the 
Einstein equations are a set of highly nonlinear equations which may exhibit dy- 
namically very complicated solutions for the space-time metric. However, in general 
relativity it is not easy to construct indicators of chaos which are gauge-invariant. 
Therefore it is reasonable to start by investigating - at first - the possibility of a 
gauge-invariant description of local instability. 

In this paper we examine an approach which aims at describing the dynamics in 
purely geometrical terms. The dynamics is formulated as a geodesic flow through 
the Maupertuis principle and a criterion for local instability of the trajectories may 
be set up in terms of curvature invariants (e.g. the Ricci scalar) of the manifold 
on which geodesic flow is generated. We discuss the relation of such a criterion for 
local instability (negativity of the Ricci scalar) to a more standard criterion for local 
instability and we emphasize that no inferences can be made about global chaotic 
behavior from such local criteria. 

We demonstrate that the Maupertuis principle implemented in the case of the 
Hamiltonian (ADM) formulation of General Relativity is equivalent to the construc- 
tion of Wheeler's superspace on which the development of three-metrics occurs along 
geodesies. Obstructions for using the superspace metric as a distance measure on 
the space of three- metrics are pointed out. The discussion is illustrated in the case 
of a particular toy-model of a gravitational collapse. 
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1 Introduction 



Due to the strong nonlinearity of Einstein's equations one is tempted to raise the issue 
of chaos in general relativistic context especially in situations where this nonlinearity is 
particularly pronounced - for example near a generic gravitational collapse to a space- 
time singularity. The property of chaos is usually understood in terms of (exponential) 
divergence in "time" of "nearby orbits" starting from slightly modified initial conditions. 
However one faces in the context of general relativity some fundamental problems which 
make the straightforward application of the machinery developed in the theory of non- 
linear dynamical systems difficult. At first, general relativity, i.e. Einstein's theory of 
gravitation, is not even a dynamical theory in the usual sense. It does not from the very 
beginning provide us with a set of parameters (describing the gravitational degrees of 
freedom) evolving in "time". Time looses here its absolute meaning as opposed to the 
classical dynamical theories where the Newtonian time is taken for granted and upon 
which all our experience in dealing with chaos is based. The division between space 
and time in general relativity comes through foliating the space-time manifold Ai into 
spacelike hypersurfaces Sf. The metric g^u on A4 induces a metric gij on St (the first 
fundamental form of S^) and can be parametrized in the form 



NiN'-N^ Nj 

Ni Qij 



where and Ni are called lapse function and shift vector respectively. Only after splitting 
the space-time into space and time (the 3-M ADM splitting) [0, H we yield the possibil- 
ity of mapping the Einstein equations to a dynamical system which resembles ordinary 
Hamiltonian dynamical systems. In general this 3 + 1 split is quite an arbitrary procedure 
and there is a priori no preferred time coordinate. 

The question we shall address in this paper is in particular the following: 

What is meant by a "nearby orbit" in general relativity? 

Let us suppose that we have performed a 3 + 1 split of the spacetime metric into a 
three-metric "(^^(yi" evolving in a time coordinate "t". What is then a natural distance 
measure | I'^^-'fi' — | between two three-metrics '^'^''g and g* at some given moment t of 
time? In the case of non-relativistic Hamiltonian dynamical systems, a standard indicator 
of "chaos" such like a principal Lyapunov exponent (cf., e.g., discussion in [|, ^ ^) is 
calculated by using an Euclidean distance measure which is naturally induced from the 
structure of the kinetic energy term appearing in the non-relativistic Hamiltonian. In 
general relativity it is, however, not obvious why we should use an Euclidean distance 
measure to assign a distance between two space-time metrics. 
In the class of spatially homogeneous spacetime metrics Q 

ds^ = -N^dt^ + -fij{t)uj'{x)u;^{x) (1.1) 

calculations of principal (maximal) Lyapunov exponents have previously been reported 



10, 12, 13| for the dynamically interesting case of the vacuum mixmaster gravitational 



^For a review of spatially homogeneous space-time metrics, see e.g. || 
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collapse. In fact, the mixmaster toy- model collapse is a useful laboratory to test ideas 
about how to characterize chaos in general relativity (see also e.g. discussions in Rugh 
r^). In the case of the mixmaster metric a;*(x) in ( |1 . 1[ ) denote the SU{2) invariant 



one-forms of the mixmaster space |1| and the three-metric is parametrized by three scale 
factors a = e", b = e^, c = , 

-fi,{t) = diag{a\t), b\t), c\t)) = diag{e^''{t), e^^{t), e'^t)) . (1.2) 

The Einstein equations for the mixmaster model are given by three coupled non-linear 
second order differential equations supplemented by a first integral constraint (cf. e.g. 
Landau and Lifshitz 0). It is tempting to treat the mixmaster collapse like any other 
dynamical system (with few degrees of freedom), integrate (numerically) the 6-dimensional 
state vector^ (a, P, 7, ar, Pr, 7r)(T) and extract "Lyapunov exponents"^ in a standard way. 
In all previous studies [1^, [T^, |I^ Lyapunov exponents were calculated by used an 
Euclidean distance measure (between two "nearby" solutions) of the formQ 

\\S(^)g,^\\' = (SaY + {5pf + {5^f (L3) 

Such choice of a distance measure is directly inspired from the study of ordinary (non- 
relativistic) Hamiltonian systems. In general relativity there is however a priori no reason 
why one should use such a distance measure to give out the distance between ^^^^f and the 
"nearby" three-metric ^^^g + 5^^^g. This distance measure has been put in artificially in 
hand and is not supported by the structure of general relativity. A distance measure which 
is induced by the structure of the (ADM) Hamiltonian formulation of general relativity is 

I \S'^'^g^J\\' = ds^ = GABSg^Sg'' = G.^MSg'^Sg"' (L4) 

where Gab = Gijki is the so-called "(mini)superspace" metric. We will here critically 
examine the idea of using the "superspace" metric Gab (or a conformally rescaled super- 
metric) as a distance measure between two three-metrics. In particular we discuss the 
idea of extracting an invariant from it like the Ricci scalar in order to give invariant (local) 
statements about the divergence of "nearby" orbits. Thus, in this paper we examine and 
extend in various ways the discussion of an approach ll^ - [^] which has proposed a 



criterion for instability of the mixmaster gravitational collapse. Szydlowski and Lapeta 
1^, Szydlowski and Biesiada proposed to look at the Ricci scalar of the manifold 



on which the mixmaster dynamics generates a geodesic flow and investigate if one could 
extract (coordinate invariant) information about instability properties and chaos of the 
gravitational collapse in this way. 

We show in sec. 2. 2 that for the case of non-relativistic Hamiltonian dynamics the in- 
stability criterion of negativity of the Ricci scalar (of the manifold on which the dynamics 



^Here t = J dt/abc denotes the standard time variable in Landau and Lifshitz and a subscript r 
means differentiation with respect to the time variable r. 

^We note that the spectrum of Lyapunov exponents is one of the most well known (and often calcu- 
lated) standard indicators of chaos in a dynamical system. 

^ Preferably, one measures the Euclidean distance constructed from the entire state vector 
(a, /3, 7, a', /?', 7') in phase space. What is equivalent - Lyapunov exponents are calculated from the 
Jacobian matrix corresponding to the state vector (a, /3, 7, a', /?', 7') and integrating up along the orbit. 
Cf. e.g. HobiU (l0| or appendix A in S.E.Rugh pi. 
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generates a geodesic flow) is related to a more standard local instability criterion, which 
examines local instability in terms of the eigenvalues of the Jacobian of the Hamiltonian 
flow. We will stress the nontrivial interplay between local and global instability criteria 
(since several authors draw false conclusions from local instability criteria). In particu- 
lar, we do not believe in expressions directly relating the average negative curvature to 



Lyapunov exponents (see, e.g., Szydlowski and Lapeta and Szydlowski and Biesiada 



We note in sec. 3 that application of the Maupertuis principle to the Hamiltonian 
formulation of the mixmaster dynamics - and the associated induction of a natural distance 
measure on the three-metrics - is exactly to implement the dream by Wheeler, DeWitt 
and others (described e.g. in pSf for the mixmaster metric) to have a superspace metric 
"Gjjfc;" (Wheeler's superspace) - with respect to which the three-metrics move along 
geodesies - and use this as a natural distance measure on the space of mixmaster three- 
geometries. In a sense such a superspace metric would be the most natural to impose on 
the configuration space of three-geometries (and appear more natural than introducing a 
completely arbitrary Euclidean metric on the space of three-geometries - as has been done 
previously in calculating Lyapunov exponents for the mixmaster collapse dynamics). 

Thus, the Ricci scalar calculated in Szydlowski and Biesiada from the Hamiltonian 



introduced by Bogoyavlenskii |^ is exactly to calculate the Ricci scalar of Wheeler's 
superspace metric and the criterion R < imply that we have local instability at least in 
one direction when we look at the geodesic deviation equation for two nearby mixmaster 
metrics investigated with Wheeler's superspace metric as distance measure. 

An obstacle to this approach is the accompanying introduction of a host of singulari- 
ties of the superspace metric which makes the dream of Wheeler troublesome to achieve 
(see also discussion in CW. Misner [^) - even in the restricted class of mixmaster three- 
geometries. Invariants calculated from the superspace metric, e.g. the Ricci scalar, inherit 
the singularities of the Wheeler superspace metric. In particular, the Ricci scalar calcu- 



lated in Szydlowski and Biesiada ||20[ has such singularities. This obstacle has also most 



recently been emphasized by Burd and Tavakol |2^. However, as a local instability cri- 
terion, R < indicates - at nonsingular points where it is defined - the local exponential 
instability in some directions of the configuration space. 

As a major obstacle, we emphasize that the distance measure induced by using the 
Wheeler superspace metric is not positive definite. The situation is somewhat similar to 
the indefiniteness of the line element in special relativity (zero distances along the light- 
cone). But in the case of Wheeler's superspace metric in general relativity there is nothing 
which prevents a trajectory (a three-metric) from crossing the null-surface. Thus, one may 
contemplate situations where a "zero" distance (between two three-metrics which need 
not be identical) evolves into a finite (positive or negative) distance. This corresponds 
formally to a (local) "Lyapunov exponent" which is "cxd" . 



The investigation is to be viewed in the light of the following question: Can one assign 
an invariant meaning to "chaos" in the general relativistic context? 

There is a problem of transferring standard indicators of chaos, e.g. the spectrum 
of Lyapunov exponents, to the general relativistic context, since they are highly gauge 
dependent objects. This fact was pointed to and emphasized in Rugh |jl2|, - and was 
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also discussed in Pullin |T5|. (See also more recent discussions in Rugh p!7|). 

One should try to develop indicators of chaos which capture chaotic properties of the 
gravitational field ("metric chaos") in a way which is invariant under spacetime diffeo- 
morphisms - or prove that this can not be done! It may very well be that there is a 
"no go" theorem, for the most ambitious task of constructing a generalization of a Lya- 
punov exponent (extracted from the continuous evolution equations) which is meaningful 
and invariant under the full class of spacetime diffeomorphisms (H.B. Nielsen and S.E. 
Rugh). One should therefore, at first, try to modify the intentions and construct instabil- 
ity measures which are invariant under a smaller class of diffeomorphisms, e.g. coordinate 
transformations on the spacelike hypersurfaces which do not involve transformations 
of the time coordinate. 

"Per time" indicators - like a Lyapunov exponent - have only a little chance of remain- 
ing invariant under diffeomorphisms. Even if we consider the restricted class of Lorentz 
transformations, i.e. space-time transformations within the special theory of relativity, it 
appears, that the time evolution of a dynamical system is slowed down if we observe it 
from a frame of reference S', which moves (is boosted) with some velocity v relative to 
the frame of reference S in which the dynamical system is at rest. Hence, the Lorentz 
7-factor will inevitably appear in the Lyapunov exponent (seen from the moving frame 
of reference S'). One notes, however, that we are able to discriminate qualitatively, in 
any frame of reference, whether or not the dynamical system is chaotic: The requirement 
that a Lyapunov exponent is greater than zero will be valid in any frame of reference 
if it is valid in one. That is, "A > 0" is a statement which is invariant under Lorentz- 
transformations - although the numerical value of the Lyapunov exponent A is not. In 
general relativity, however, "A > 0" is not an invariant statement under diffeomorphisms 
more general than Lorentz-transformations [|T^, [T^, [T^, |TB| . 



A route of progress may lie in pointing to indicators of purely geometrical nature 
("fractal dimensions", etc). We shall here explore into a picture which is purely geomet- 
rical and stems from an old construction known since times of Maupertuis and Jacobi — 
namely that in which Hamiltonian dynamical systems are mapped into geodesic flows. In 
this approach the entire complicated dynamics comes out from purely geometrical prop- 
erties of a single "object" — the manifold on which the geodesic flow is generated. One 
may consequently dream that also such basic properties like integrability or chaos may be 
somehow captured in purely geometrical (and hence timeless and invariant) terms. For 
example one may hope to describe local instability of trajectories in terms of curvature 
invariants (e.g. the Ricci scalar) since it is curvature that determines the behavior of 
close geodesies (through the geodesic deviation equation). We shall however also point 
out several obstructions one has to face when implementing such a line of thinking to 
general relativistic problems. 



2 The dynamics as geodesic motion by application of 
Maupertuis principle 

It is indeed a remarkable fact that a large class of Hamiltonian systems can be presented 
in such a way that Hamilton's equations acquire formal resemblance to geodesic equations 
in Riemannian geometry. 



5 



The phase space trajectories x{t) = (p(t), q{t)) with respect to time parameter t and 
corresponding to the Hamiltonian 



H = H{p, q) = Y: \a'^V^VJ + V{q) (2.1) 

may be mapped into a geodesic motion on the configuration space manifold. It is accom- 
phshed in a way which displays a complete analogy with the mapping of the evolution of 
three-metrics into that of geodesic motion in Wheeler's superspace (see next section). 
The equations of motion for the q coordinates may be written as 

q^+T\,q^q' = -a^'^-^ (2.2) 

where T\g is the Christoffel symbol calculated with respect to a^j metric. Due to the force 
term this is, obviously, not a geodesic equation. It is simply the Newton's second law 
restated. The momentum variables are just linear combinations of velocities pi = aijqK 
Transformation to a geodesic motion (i.e. free motion in a curved space) is accomplished 
in two steps: (1) conformal transformation of the metric Oj^, and (2) change of the time 
parameter along the orbit. In fact only the first step is crucial the second one merely intro- 
duces the affine parameter along the geodesies. More explicitly we equip the configuration 
space with the (super) metric 

Qi, = 2{E - Viq))aij (2.3) 

(note that a^j is read off from the kinetic energy term in the Hamiltonian and (in general) 
is allowed to vary as a function of the configuration space variable, aij = aij{q)) ). 
With respect to the metric ( ^.3] ) and the time parameter t it is not easy to see that the 
orbits are geodesies since there is a term appearing on the right hand side of the equation. 



where T^jf, now denote the Christoffel symbols associated with the (super)metric 
However, if we (re)parametrize the orbit = q^{s) in terms of the parameter s defined 
via 

s = J2{E-V)dt (2.5) 

the orbits will become affinely parametrized geodesicsQ, i.e. the configuration space vari- 
ables q = (g*) satisfy the well known geodesic equation 



^Let us recall from differential geometry that the most general form of the geodesic equation (when the 
geodesic is parametrized by v parameter) reads q * + F^j.g-'g = h{v)q ' (prime denotes differentiation 
with respect to v) where h{v) is an arbitrary function of v. Then by solving the differential equation 
^ + h{v)^ — one can construct a new parameter — the so called affine parameter s such that 
9* + Tj-fe?"''?'^ = (overdot here denotes differentiation with respect to s). It is evident that this new 
parameter is defined up to linear transformations s — > As + B (A, B are constants) what justifies its 
name. 
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with no force term on the right hand side. 

There is no "magic" in this, it is merely a rewriting of the Hamiltonian equations. 
The information about the original force acting on the particle (as described by the 
potential V{q) in the Hamiltonian has been encoded entirely in the definition of 

the (super)metric ( ^.31 ) and the definition of the new parameter s in ( |2.5| ) parametrizing 
the orbit. The entire procedure involves the following principal steps: 

1. At first a metric a^j is read off from structure of kinetic term in the Hamil- 
tonian. 

2. We observe that the trajectory {qit)} is not a geodesic with respect to this 
metric. 

3. Transformation of the non-geodesic motion to affinely parametrized geodesic 
motion is obtained via two steps: 

(1) Conformal transformation of the metric aij — > 2{E — V{q))aij 

(2) Rescaling of the time parameter t ^ s = J 2{E — V{q))dt 

It is, however, not often that this description is used in the discussion of Hamiltonian 
dynamical systems, for example in the discussion of instability properties and chaotic 
behavior of such systems. We contemplate several reasons why this approach is so rarely 
used: 

1 . The metric Qij in ( |2.3| ) and the parameter s in (|2.5| ) are singular in "turning" 
points where the kinetic energy term J2 \'^^^ ViVi is zero and we thus have 

E - V{q) = . 

Therefore we have to consider the motion away from the "turning" points. 
This may be a more severe restriction for systems with few degrees of free- 
dom, e.g. the harmonic oscillator with a single degree of freedom, than for 
systems with many degrees of freedom, e.g. for N self-gravitating bodies in a 
gravitational field, to which the virial theorem applies, and thus J2 ^O'^'^PiPj = 
is rather unlikely. 

2. The procedure of mapping to geodesic motion (via the Maupertuis prin- 
ciple) is applicable to any time independent Hamiltonian system of the form 
(|2.1| ) where the potential V{q) is not allowed to depend on momenta p. I.e. 
the application of Maupertuis principle is not straightforward if we for ex- 
ample have a mechanical system coupled to an electromagnetic field with the 
Hamiltonian 

H=^ip--Aiq)r + Viq). 
2m c 

®Note that the above picture comes quite naturally from the Maupertuis- Jacobi least action principle 



6S = 6 / ^E~V{q)Jaij dq'dqi = 

Jq' 

and its formal resemblance to the variational formulation of geodesies in Riemannian geometry as curves 
extremalizing the distance. We have sketched the step-by-step derivation of the geodesic motion in order 
to see the step-by-step analogy with the derivation of three-metrics as geodesies in Wheeler's superspace 
in the next section. 
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The same objection holds for mechanical problems in noninertial frames, for 
example the problem of motion of a star in the field of a rotating galaxy EH . 



In spite of these principal limitations of the Maupertuis approach, it has been explored 
in various contexts. In his seminal work |^ on foundations of statistical mechanics 



Krylov has suggested that viewing the dynamics of an N-body system as a geodesic 
flow on the appropriate manifold may provide a universal (applicable to a very large 
class of interaction potentials) tool for discussing relaxation processes. Inspired by works 
of Hedlund and Hopf [0 Krylov conjectured that the average separation of trajectories 
evolves according to the sign of the Ricci scalar R on the configuration space accessible for 
the system. An implementation of this viewpoint have also been attempted in the context 
of general relativity by Szydlowski, Biesiada et al. ITS, ITP|, |2D|. As we shall demonstrate. 



Wheeler's original dream to construct a superspace in which the three metrics move along 
geodesies is exactly to implement the old idea of Maupertuis principle in general relativistic 
context. We shall however first - in the next two subsections - comment on the status of 
the R < criterion and its relation to more standard criteria for local instability. 

2.1 The Ricci scalar R for the "super" metric and < as a 
(strong) local instability criterion. 

Contemplating the local instability properties, i.e. how nearby orbits behave, it is nat- 
ural to consider the geodesic deviation equation which describes the behavior of nearby 



geodesies (|2.6| ) 



This can be derived in the usual manner by subtracting the equations for the geodesies 
and respectively or simply by disturbing the fiducial trajectory {pi(t), q^t)), 

Piit) = Piit) + r]i{t), 

q\t) = q\t)+C{t) (2.7) 

and substituting this directly into Hamilton's equations. In this way we also arrive mo- 
mentarily, though tediously, at the geodesic deviation equation for the configuration space 
variable ^ (the separation vector). 



Ds^ 



(2.8) 



Here u = Dq/Ds is the tangent vector to the geodesic, $, is the separation vector orthogo- 
nal to u. Note, that the covariant derivative D/Ds and Christoffel symbols are calculated 
with respect to the (super)metric (|2.3|). This (super)metric induces the natural distance 
measure 

\m'=9^Je^' =2{E-V{q))a,je^^ 

on the configuration space. From the geodesic deviation equation ( ^.81 ) one arrives at the 
inequality (see e.g. p^ ) 

^iiiir>-2i^.diiir (2.9) 
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where the K^-^ denotes the sectional curvature in the two-direction u A ^ . 

Noting that the Ricci scalar may be expressed as a sum of the principal sectional curva- 

tures,[] 

R — Tr K — ^ ( -^^^"''^^^"^^ Sectional \ ,^ 
^ \ Curvatures 1 

the inequality (|2.9| ) especially implies that if 

R<0 (2.11) 

then a principal direction ^ exists along which we have local exponential instability. 

By first calculating the Christoffel symbols T^j from the metric ( |2.3| ) and subsequently 
substituting this into the standard formula for the Ricci scalar R one obtains the following 
general expression, 

' - ^ £ - w"" - ^li-} 

corresponding to a Hamiltonian system with n degrees of freedom. 

N.S. Krylov was - as far as we know - the first to emphasize the use of i? < 
as an instability criterion in the context of an body system (a gas) interacting via 
Van der Waals forces, with the ultimate hope to understand the relaxation processes in a 
gas. In his toy-model study he found that indeed i? < in the accessible domain of the 
configuration space and concluded that this provides an explanation for the relaxation as 
a consequence of dynamical mixing.^ Another investigation in the same spirit, have been 
performed by Gurzadyan and Savvidy |^| for the case of collisionless stars interacting 



via gravitational forces, also with the purpose of understanding the relaxation properties 
of such a stellar system. 

The system of collisionless stars which was investigated by V. Gurzadyan and G.K. 
Savvidy |2^ and Kandrup |Q had R < whenever the number of stars in the system 



exceeded two but not on a domain which was compact. Kandrup P3, however, argues 



''The principal sectional curvatures are defined, as is well known, as the extrema of the sectional 
curvature K^-^ as function of the wedge product u A ^. 

*Note that if there is one direction ^ (for a given u) in which there is exponential instability then 
there is also a continuum of nearby exponentially unstable trajectories. Due to smoothness of Kui^ as a 
function of two-direction the sectional curvature for nearby directions ^ is also negative. 

^The investigations by Krylov was strongly influenced by seminal results of Hedlund and Hopf 
where the geodesic motion on the Lobachevsky space was proved to be mixing. Whereas the approach 
by Krylov is intuitively appealing it is however not allowed to transfer the results obtained for the 
Lobachevsky space to the case where sectional curvatures are negative only in average and not at every 
point and in every two-direction u A ^. Moreover, Krylov disregarded the problem of compactness of the 
configuration space manifold which is important for making inferences about mixing. 

^"Physical motivation of this study comes from the observation that stellar systems (globular clusters 
and galaxies) are apparently in a well relaxed state that is reflected in regularity of their shapes, velocity 
dispersions, surface luminocities etc. On the other hand the most obvious relaxation process - namely via 
binary encounters, provides a relaxation time that is greater than the Hubble time [^ . As a possible way 
out of this puzzle Lynden-Bell |]3^ has proposed the so called violent relaxation mechanism. Gurzadyan 
and Savvidy attempted at looking at the relaxation process in stellar systems from Krylov's viewpoint. 
The idea was that if they could show that a self-gravitating N-body system was a K-systcm it would 
indicate that the system approached the equilibrium exponentially fast. 



9 



that the self-gravitating system is "effectively bound" i.e. the phase-space is compact at 
least for sufficiently small times. If one waits long enough the effects of "evaporation" of 
stars from the system come into play, so the phase-space is actually infinite. Let us 
emphasize that the negativity of the Ricci scalar provides only a sufficient condition for 
local instability and it is thus not possible to make any rigorous statements about chaos 
from local instability properties alone as it will be discussed in more detail in the next 
subsection. 



2.2 Remarks on the instability criterion R < and its relation 
to more "standard" instability criteria. 



According to Szebehely there are more than fifty distinct criteria of instability! Many 



of these are related to each other (though not all connections are worked out in a trans- 
parent way). Some of the criteria are even in conflict with each other. 

Let us emphasize a relation between two important local criteria, the criterion of having 
a negative Ricci scalar R = R{Maupertuis) < (corresponding to the formulation via the 
Maupertuis principle) and the more usually applied criteria connected to the eigenvalues 
of the Jacobian of a Hamiltonian flow. 

-R < is a sufficient criterion for local instability — but it is not a necessary criterion 
and very often it fails to say anything of interest about our dynamical system. For 
instance, in the case of a simple class of Hamiltonian flows 

H = E = ^{pI + pI) + ^{x' + y') + ^xV (2.13) 

it fails to mirror the transition from almost integrable to almost completely chaotic (almost 
K-system) as the energy E increases. A calculation of the Ricci scalar, according to 
formula (|2.12|) or ( p.l7| ) shows that it is always positive irrespective of the Energy E 
(since AV{q) = 2 + x"^ + > 0). Despite this, the model is almost a K-system in 
the hmit for large energies (i.e. when the nonlinear term becomes very large), see 

discussions in P. Dahlqvist and G. Russberg |^ and references therein. 



A standard criterion for local instability is to look at the eigenvalues of the Jaco- 
bian of the Hamiltonian flow. (However, wrong conclusions as regards global chaos are 
often drawn from it.) The line of thinking starts at disturbing the fiducial trajectory 



^-'^Had the manifold been compact it would - of course - still not be possible to deduce that the system is 
a K-system. The K-system property holds for geodesic flows on compact manifolds of negative curvature, 
i.e. manifolds whose sectional curvature Ku;^ is negative at every point and in every two-direction u A^. 
On the other hand it is a general property of Hamiltonian systems that sufficiently close to the physical 
boundary of motion there exist both positive and negative sectional curvatures |35| . This means, in 
particular, that mapping a Hamiltonian flow to a geodesic flow via Maupertuis principle can never be 
used to show that the system is a K-system. We note, however, that to be a K-system is a very strong, 
highly non-generic property to satisfy. If a system is "less" than a K-system, it may - of course - 
nevertheless be "highly" chaotic (and unstable almost all over in the phase space). 

^^In an attempt to predict the onset of dynamical chaos (a key ingredient of which is the exponential 
instability of adjacent trajectories) Toda |3|] and Duff and Brumer proposed the criterion based on 
evaluation of the Gaussian curvature of the potential V(q). 
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q\t) = q\t)+C{t) (2.14) 

and subsequently investigating the evolution of the disturbance vectors ^ and rj. The 
stabihty of motion is determined by integrating along the trajectory (see later) the time 
dependent matrix (the Jacobian matrix of the Hamiltonian flow): 

where the Hessian of the potential V{q) is evaluated along the reference trajectory making 
the J- matrix time dependent. The inconvenience of the time dependence can be overcome 
by replacing the time-dependent phase-point q{t) (moving along the trajectory) by a fixed 
phase-space coordinate q resulting in 

e = r/^ = a%- (2.16) 
f dq^dqi ^ 

This step is sometimes claimed "mathematically dubious" (cf. also, e.g., M. Tabor 
[PI). Anyway, the J-matrix is now time- independent and the stability of the autonomous 
system ( |2.16|) is determined from the eigenvalue problem 

det [J - Al] = 

If at least one eigenvalue has a positive real part we have exponential growth (locally) 
of the disturbance vector in one direction of the phase space. For example, in the case of 
a 2-dimensional Hamiltonian we have the two pairs of roots 

1/2 



. , ,d^v d'^v, rw ;tww tipv 





^Q2qi Q2q2 y ^2^1 ^2^2 Q2qi Q2q2 dq^dq 



In particular, if det{-^^) = ^^^^ — ( 3^13^2 )^ < (i.e. negative Gaussian curvature of 
the potential V) a pair of roots becomes real and we have local instability. 

The R < criterion implies the J-instability above|^ (whereas the converse is not 
true). In this sense the i? < is a stronger instability criterion than the J-instability. For 
example, for the two dimensional Hamiltonian system (with aij = 6ij) the Ricci scalar 
( gl2D reads: 



R{Maupertuis) = - \ /\\n{E -V{q)) 

2{h - V{q)) 



[AV{q) + ^^^[j^^l ) (2.17) 



2iE-Viq)y' E-V{q) 



where A = d /d{q^) is the Laplacian (if a^j = 5, 
Thus from < follows that /\V < which in turn implies the existence of eigenvalues A 



^^This is not surprising from simple considerations: If we have instability along some direction in the 
configuration space (as captured by the R < criterion) then we also have instability at least in one 
direction in phase space. 
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with a positive real part irrespective of the sign of the Gaussian curvature of the potential. 
On the other hand the J-instability criterion (sometimes called the Toda-Brumer criterion) 
may capture the change of local instability properties of the Hamiltonian flow even in cases 
when AV > and hence R > . For example, the instability criterion R < fails to 
capture the transition from KAM-integrability to chaos (as the energy E increases) in 
the class of Hamiltonian flows ( p.l3| ), whereas the study of the Jacobian of the flow (the 
Toda-Brumer-Duff criterion) gives a threshold value E = E^. (onset of "local instability" ) 
for the energy. Thus, a small calculation will verify that for energies 

E> Ec = 3/2 

there are eigenvalues of the Jacobian J{x) with a positive real part. Let us emphasize, 
what is well known, that it is not possible to make conclusions about chaos from this local 
criterion. Global instability (chaos) is defined by looking at the spectrum of Lyapunov 
exponents: Consider the set of first order differential equations, x = f{x), x = (p, q) G 
i?" (the Hamiltonian flow) . The Jacobian J of the flow / is the n x matrix J = {df / dx) . 
In order to define the Lyapunov spectrum of characteristic exponents connected to a given 
trajectory x(t) with initial conditions x{0) one first define the stability matrix, 

^ = P-i«> 

which satisfies the differential equation 

Mij = J{x{t))ikMkj , M(0) = 1 (2.19) 

The solution to this differential equation is called the exponential of the Jacobian under 
the flow, symbolically denoted by 

M{t) = f exp( /* J{x{t))dt)M{0) (2.20) 
Jo 

where T is the time ordering operator (the J's do not commute, being n x n matrices) .[^ 
In order to determine whether a trajectory is stable (regular) or not, one investigates 
the growth of the stability matrix M = M{x{0),t) = M{t) with time. More explicitly, 
consider A{t) = {M\t)M{t)y/'^ (cf., e.g., Eckmann and Ruelle, § p. 630). Let the n 
eigenvalues of the n x n matrix A{t) (which are real and positive) be denoted di{t) > 
d2{t) > ... > dn{t). The Lyapunov functions are defined as Xj{t) = t~^ \ogdj(t) and the 
associated spectrum of Lyapunov exponents 0, 

A, = hm X,{t) = hm |i logrf,(t) j . (2.21) 



T'exp(Jg J{x[t))dt) may be calculated as TliniAt^o !!(-'- + •^(^(^)^^)-^(O) the product being taken 
along the integrated trajectory {x{t)} in time steps At 0. 

^^Theoretical criteria for convergence of this limit, which holds under rather general assumptions of 
the flow /, are discussed in e.g. Eckmann and Ruelle |Q. The maximal characteristic exponent \max is 
the easiest to find. But this is also sufficient to detect the exponential divergence of nearby orbits, called 
"sensitive dependence on initial conditions" . Per definition, a dynamical system displays a stochastic 
chaotic behavior in a region of the phase space i/its maximal Lyapunov characteristic exponent is positive 
for trajectories in this region, cf. Benettin et al. 



12 



In the very special case where the Jacobian is constant in time, J{x{t)) = Jq we get 
M(t) = exp( Jq 't)- If A is an eigenvalue of Jq we have that e^* is an eigenvalue of M{t). 
It follows, that A is one of the Lyapunov exponents and local J-instability is equivalent 
with global instability (Lyapunov exponents). 

In general, however, no implication exists between local and global instabilities. This 
is illustrated by a simple model example. Consider, e.g., a two-dimensional map 

X {ABABABAB...)x , x e 

where a 2 x 2 matrix A is multiplied on the vector every second time, and the 2x2 matrix 
B every second time. 

If we choose A = ^ and B = ( \ ^ \ \ E R, both are locally stable, since the 



,0 IJ \X 1. 

eigenvalues are equal to 1. AB has eigenvalues 1 ± A and globally the system is unstable. 

On the other hand if we choose A = ( ?^ and B = A~^ = (^^J^ V where 



a; ^^^^ = ^ = I l/A, 
A > 1 and < A < 1, both A and B gives local instability while globally the system is 
stable {AB = 1). 

The same non-trivial relationship between local and global instability criteria holds for 
continuous Hamiltonian flows: Orbits 7 may build up a "global instability" (as captured 
by the exponentiated Jacobian along the orbit, cf. formula (|2.2CI|) ) despite they are stable 
everywhere in the phase space with respect to the local instability criterion. For example, 
if we consider the simple Hamiltonian toy- model ( p. 131) , for energies E < Ec = 3/2 there 
are still many orbits which are globally unstable, despite the fact that there are no positive, 
real eigenvalues of the Jacobian in any accessible phase space point for the trajectory. 



Moreover, we may have (globally) stable orbits (Dahlqvist and Russberg ||3^) which 
are everywhere locally unstable. (Consider the Hamiltonian ( |2.13| ) in the limit E 00. 
The Jacobian has real, positive eigenvalues in any point of the accessible phase space. 
Nevertheless, there is an orbit of period 11 which is globally stable, cf. Dahlqvist and 



Russberg |37 



Thus there is a non trivial interplay between the local instability (as mirrored in the 
eigenvalues of the stability matrix (|2.18|) ) and the global, exponentiated Jacobian ( |2.2CI| ), 



yielding the "Lyapunov exponents" along the orbit. These remarks are well known. We 
repeat them here to comment on what one may conclude (or not conclude) from local 
criteria of instability like 

R{Maupertuis) < 0. 

There has been some confusion concerning this point in the past. 

Especially, we can not from R < (local instability) deduce that a Lyapunov exponent 
is positive (global instability). 



As already noticed in the introduction there were attempts [19, EQ, 21] to relate the 



R < criterion with Lyapunov exponents. In particular it was argued that the quantity 



A = hm^^^J ^ - r2{E - V{q{t)))dt (2.22) 
y n[n — 1) r Jo 

played a role of the principal Lyapunov exponent. The Ricci scalar R is meant here as an 
average measure of local divergence of trajectories. The motivation of such claims (see 
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[ p!9| , pO| , pT| ) had its roots in the following facts: If we take the geodesic deviation equation 
(^]8|) and rewrite it as a gradient system (see appendix 1 in Arnold fS^), i.e. 

where the "potential" ytt(^) has been expressed by sectional curvature: Vu{^) = -Ku-^g{^,^). 

Then by averaging the geodesic deviation equation over the orientation of the geodesic 
(i.e. at a given point we take all possible geodesies passing through this point, average 



over the bivector n A ^ and denote this average by < >) one arrives at the formula 



Ds"^ nin — r 



<e>, 



from which one may argue for the quantity (p.22|) bearing in mind that the s-parameter 



is related to the t-time by (|2.5| ). The reasoning is flawed, however, in several aspects. 
One cannot replace the integration of a full Jacobian matrix, cf. equation ( p.20| ), with the 



integration of a Ricci-scalar which is just some average quantity (trace over the sectional 
curvatures). The negativity of the Ricci scalar, i? < 0, which is an initial assumption 
for ( |2.22| ) to apply is only a sufficient criterion of local instability. Hence in the case of 
the Hamiltonian ( p^.l3 ), for example, we have i? > whereas the Lyapunov exponent 



A > for most trajectories. Then apart from the fact that local instability cannot be 
simply translated into global chaos also the meaning of "averaged deviation equation" 
is questionable. Even if the Ricci scalar R is traced (numerically) along trajectories one 
may give examples of systems with positive Kolmogorov entropy (understood, here, as 
the sum over positive Lyapunov exponents) for which the quantity ( p.l2| ) is either positive 



or negative (see e.g. discussion in [^). The above critique apply to some earlier papers 
[0, |2y, |21[] and also to Burd and Tavakol who adopted the expression ( |2.22| ) in p7[] . 

We would also like to question the validity of the formula (38) in Gurzadyan and 
Savvidy in which the Ricci scalar is related to a relaxation time of the dynamical 



system consisting of collisionless stars interacting via gravitational forces. 

In order not to close this section in too pessimistic mood let us stress that the 
Maupertuis-Jacobi procedure of reducing the dynamics to geodesic flows is a powerful 
tool for geometrizing the dynamical systems[3 (provided one is bearing in mind its limi- 
tations, see also p3[] ). 

A serious problem in general relativity is to define a metric space in which our space- 
time metric (which is itself a metric space) is a phase space point and where the divergence 
of trajectories makes sense. The next section shall reveal the intimate connection between 
such metric structure called the superspace and the Maupertuis-Jacobi principle imple- 
mented to the General Relativity in its Hamiltonian (ADM) formulation. 

The geometrical setting arising from the Maupertuis-Jacobi principle makes this picture suited for 
investigating integrability of Einstein's equations as a consequence of existence of Killing vectors and 
tensors in the minisuperspace |Q. Also, the fact that the Maupertuis principle involves variations at a 
fixed energy makes this formulation useful for the microcanonical ensemble and proves useful in deriving 
thermodynamical properties of gravitating systems as advocated by Brown and York [Q . These possible 
applications of the Maupertuis principle are beyond the scope of the present paper and the interested 
reader is referred to a review paper |23] and the references therein. 
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3 Wheeler's superspace: Maupertuis principle im- 
plemented in the context of general relativity. 

An original dream by Wheeler, DeWitt and others is to consider the dynamics of the three 
metrics ^^^g as geodesies on some manifold called superspaceQ equipped with the metric 
tensor Gijki which is named the supermetric. This supermetric Gijki induces a norm on 
the space of three metrics ^^^g, 

WSgabW = J d'x^ G'^'''6g,,6gki (3.1) 

and one may measure the distance between two three-metrics ^^''g and ^^^g with respect to 
the Gijki tensor. (In practice it may be difficult to find such global distances. Like in the 
case of finding distances between Copenhagen and Warszawa, say, one has to minimize 
over all paths connecting the given two points). 

What can one choose for the distance measure Gijki ? At first there are many possible 
choices of metrics on the space of three metrics. For example (cf. DeWitt |^) one could 
choose a metric like 

which has the property that the distance between two three-metrics is zero if and only 
if the three-metrics are identical. However, another distance measure (supermetric) is 
suggested from the ADM Hamiltonian formulation of general relativity. It is read off 
from the kinetic term in the Hamiltonian constraint, 

^■<-\{<f-9'^'^R = ^ (3.3) 



(cf. e.g. Misner and references therein) associated with the constrained Hamilto- 
nian formalism which may be set up for the Einstein equations (see e.g. discussion and 
references in Teitelboim ^^)- More explicitly. 



1 



T^k? = -T^ig'^g^^T^ijT^ki + g'^g^'^T^ijTrki - g'^g'^^rtijiTki) = G'^'^^TTijirki 



2' 2 
and we thus arrive at 

^ ^^gikgjl ^ gUg^k _ g^ ^ ^ ^^ ^^^ 

Note that the supermetric, in this notation differs by a conformal factor (y^) from 



DeWitt's expression This is allowed since the ADM action is invariant with respect 
to conformal transformations. 

Let us now see how we may arrive at the geodesic equation for the three-metric '^^'>g. 
It turns out that it can be done in a way which is completely analogous, step-by-step, 
to the previous section dealing with the similar question for an ordinary (nonrelativistic) 
Hamiltonian system H = ^aijp^'p^ + V{q). 

The procedure involves the following principal steps: 



^^Abstractly the "superspace" is meant to be the space of all three-metrics modulo difFeomorphisms 
on the three-space, i.e. the quotient space S{A4) = Riem{A4)/ Dif f{M). 
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1. The metric Gab is read off from tlie ADM Hamiltonian. 

2. Observe tliat we do not have geodesic equation for {^g}, with respect to 
metric Gab- 

3. Transformation to (affinely parametrized) geodesic motion of {^g} is ob- 
tained in two steps: 

(1) Conformal transformation of the metric to Gab = T^Gab- 

(2) Rescahng of the A parameter, A — > A. 

1. As we have aheady noted, from the structure of the Hamiltonian 

H = lG^^Mkl)^''^'' - lg^''>R = \gab7^^tt^ - \g '^'^R = O (3.5) 

we read off the first candidate for a metric on the configuration space (the space of three- 
metrics) Q 

Gab = G(.ij)(^ki) = -{gikgji + gugjk — '^gijgki)- (3.6) 

This step corresponds to reading off the "metric" a*-' from the Hamiltonian 

H = ^a^^p.pj + V{q) = E 

in the previous section. 

2. We observe that we do not have a geodesic equation for the trajectory '•'^^(7(A) 
There is a "force term" on the right hand side: 

<Pg^ . dg^ dg^ _ 1 dU _ 1 ^ABdi/^ 

d\^ +^^^dA d\ "2^ dg^~2^ dgB ^^'^^ 

where 

-pA ^r^AD/ 9GBD , dGpc OGbc^ /o o\ 



Equation ( |3.7| ) may be one-to-one translated to the similar expression (|2.2| ) from sec. 2, 
where we also have a force term, —a^^dV{q)/dq'^, in the form of the metric (read off from 
the kinetic term) multiplied by the gradient of the potential. The role of the potential in 
this general relativistic context is played by the quantity 

V = -lg('^R^-^n . (3.9) 

In analogy with section 2 (Maupertuis principle) we want to make the three-metrics move 
along geodesies ! 



"'^^The covariant form of the supermetric Gijki actually differs from its contravariant counterpart. This 
is the consequence of a natural duality of these two tensors: 



"'^^By A we denote the parameter which parametrizes the evolution of three-metrics — the so-called 
supertime. For more detailed discussion see C.W.Misner |25 
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3. Transformation of the non-geodesic motion (above) of ^^^g to geodesic motion is 
obtained in two steps: 



1^1) First we make a conformal transformation of the metric 



GAB = nGAB , G^^ = iG^^ (3.10) 



This conformal transformation of the metric (|3.10|) is analogous to equation (|2.3|) in section 
2, i.e. the conformally rescaled metric, 

g,, = 2{E - V{q))ai, 

since 

Gab = TZGab = i-2V)GAB = 2{E - V)Gab (3.11) 

(note that the energy E is zero in the general relativistic context by virtue of the Hamil- 
tonian constraint). 

Analogously to sec. 2 conformal rescaling of the metric alone is not sufficient to map the 
evolution of three-metrics into affinely parametrized geodesies. 

(2) However, by rescaling of the A parameter, 

d\ = 2{E - V)dX = -2VdX = TZdX = g'-^^R dX (3.12) 

we obtain now the important result (cf. also e.g. Misner p5|) that with respect to this 
new parameter A and the conformally rescaled metric Gab the three-metric ^^''g is now 
an affinely parametrized geodesic (in this "Superspace" ) 

-^ + Tj,r^^ = ^. 3.13 
dX' dX dX ^ ' 

Like in the previous section there is no "magic" involved (we have "magic without 
magic"). Information about the potential ^^^R has simply been encoded completely in 
the mathematical definition of the Superspace-metric G with respect to which the evolu- 
tion of '^^'>g then becomes geodesic motion (affinely parametrized if one uses the rescaled 
parameter A). 

Note, that instead of Misner's analogy with a free particle in special relativity, with 
Hamiltonian H = ^{ri^'^p^Pu + ^^) (cf. Misner [^, p. 451) we have rather emphasized 



here the complete one-to-one correspondence between construction of the Superspace and 
dynamics of a non-relativistic particle with the Hamiltonian H = ^a^^PiPj + V{q) reduced 
to geodesic flow by virtue of the Maupertuis principle. 

Whereas we would not be surprised if the analogy above is to some extent well known 
for those who understand the construction of Wheeler's superspace, we have not seen 



the analogy stated so clearly and it puts the previous investigations |jT8|, 0, ^ |2T|, ^ 
along this route in a somewhat new perspective. Namely they can be viewed as exploring 
the geometric structure of the superspace aimed at investigating local instabilities in the 
evolution of three-geometries. Dreaming of this purely geometrical picture of dynamics 
one has nevertheless to face several problems: 
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1. There may exist (and in fact it is not of rare occurrence) points where the potential 
term TZ is zero. At such points the Supermetric ( p.lO|) as well as the new parameter 
A = / IZdX are not well defined and present at first an obstruction to this beautiful idea 
of mapping the evolution of three-metrics into geodesic motion w.r.t. this Superspace. 
Of course as already mentioned in sec. 2 it is merely an artifact which stems from the 
conformal transformation ( |3.10| ) and nothing singular happens at those points to the 
dynamics. This fact has been well known already to Misner P3| . It has also recently been 



emphasized by Burd and Tavakol p7 



2. If we are to uphold the dream in the strongest version of having the three-metrics 
as geodesies in superspace it should be the conformally rescaled "super" metric Gab 
which should be called the "super metric" (since that is the one which generates geodesic 
motion of the ^g) and not the Gab = Gijki metric, which is more frequently used as a 
"superspace metric" . This imply, however, that the (mini) superspace metric would differ, 
for example, from one minisuperspace model of some Bianchi type to a minisuperspace 
model of another Bianchi type - since the conformal rescaling depends on the specific 
space-time metric considered. 

3. The Supermetric Gab induces a natural distance measure on the three-geometries 
in the relevant superpace. However, it should be noted, what is very well known, that 
Gab represents a pseudo-Riemannian geometry (with signature ( — 1-++++)), therefore 
we may have distances ds"^ which are space-like, null or time-like. (This indefiniteness of 
the distance measure is particular to the application of Maupertuis principle in general 
relativity and the natural distance measure induced by this procedure. Application of 
Maupertuis principle to non-relativistic Hamiltonian systems give rise to positive definite 
distance measures). Therefore in general relativity we may have the unfamiliar situation 
that two three-metrics which are at zero distance with respect to the Wheeler's superspace 
metric Gab naay evolve into a finite (positive or negative) distance. This corresponds 
"formally" to a (local) "Lyapunov exponent" which is "cxo" . 



In analogy to the non-relativistic case but bearing in mind all the limitations we may 
contemplate the behavior of close trajectories (mapped into geodesies) starting from the 
geodesic deviation equation: 




The idea is to use 

^ jn-i) ^ f d'n {n - 6) on an 

Ab\ dg^dg^ 4 dg^ dg^ 

as a geometrical (and hence coordinate invariant) local instability criterion for the three- 
metrics evolving according to the Einstein's equations. If i? < then there exist at least 
one direction of perturbing the three-metric g^ ^ g'^ + 5g^ so we have a local exponential 
amplification of the metric perturbation in that direction - measured with respect to the 
Wheeler's superspace metric as a distance measure on the three-metrics. 
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In the next section we shall apply this in a concrete example of the orbits of the 
mixmaster minisupcrspace gravitational collapse. 
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Supermetrics (Maupertuis principle) and local instability criteria 
in non-rel. systems and in general relativity 


A Hamiltonian for a 

non-relativistic 
mechanical system 


ADM-Hamiltonian 
in General Relativity 
(say, the mixmaster collapse) 


1. Read off tlie first candidate 
of a supermetric "aj/' from the structure 
of the kinetic term in the Hamiltonian 


1. Read off the first candidate 
of a supermetric "Gab" from the structure 
of the kinetic term in the ADM-Hamiltonian 

H = E iGAB-rr^rr^ - \g ^^^R = 


2. We observe that the trajectory {(^(t)} 
do not obey a geodesic equation of motion 
w.r.t. metric a*-': There is a force term 
-a^^dV{q)/dq^ 
on the right hand side of (^.21) . 


2. We observe that the trajectory Ifi'"^! 
do not obey a geodesic equation of motion 
w.r.t. metric G^^: There is a force term 

on the right hand side of (3.7). 


3. Transformation to geodesic motion 
is obtained in two steps: 
(1) Conformal transformation of the metric 
a,-i 2(E — V(Q))aji 
(2) Rescaling of the time parameter 
dt-^ds = 2{E - V{q))dt 
and we obtain the equation (p. 61). 


3. Transformation to geodesic motion 
is obtained in two steps: 
(1) Conformal transformation of the metric 

Gab = 2(E - V)Gab = [q ^^^R)Gab 
(2) Rescaling of the A parameter 
dX = 2{E - V)d\ = (g ^^^R)dX 
and we obtain the equation (|3.13|). 


The Ricci scalar may be calculated 
from the resulting "supermetric" 
and the criterion 
R<0 

is used as a local instability criterion 
(N-body systems etc., Krylov etc.) 


The Ricci scalar may be calculated 
from the resulting "supermetric" 
and the criterion 
R<0 

has been attempted as a 
local instability criterion 
in the context of general relativity 
(Szydlowski, Biesiada etc.) 



Table 1 : summarizing the analogy between the Maupertuis principle for non-relativistic Hamil- 
tonian systems and Wheeler's superspace (Maupertuis principle) for general relativity. 



20 



3.1 Application to the mixmaster gravitational collapse 

For the mixmaster homogeneous toy-model collapse with the three-metric, 

where 7ij(t) = diag{a'^{t),b'^{t),c'^{t)) and a,b,c are the scale factors of the metric, we 
have in the ADM-variables, 




/ ^3 ^3 
6 6 




(3.14) 



v 

the following form of the Hamiltonian 

n = \{G^%APB -n) = \{-pI + pI +pl + e-^"(y(/3+,/3_) - 1)) (3.15) 

where the mixmaster three-curvature potential reads 

n = _e-4f^(y(/3^,/3_)-l) (3.16) 

= -e-^"[-e^^+(cos/i(4V3/3_) - 1) + -6"^^+ - -e-'^+cos/i(2V3/3_)] 
3 3 3 

The (mini)superspace variables are thus = {g^,g~^,g~) = {^,13+. {A.B — 
1,2,3) and the momenta are Pa = GABdg^ /d^ S^^i^SPn = —dQ/dX ,pf3^ = df3^/d\ , pp_ ~ 
dP-JdX. We immediately read off the first candidate for a "super" metric from the kinetic 
term in the Hamiltonian 

Gab^G^'' ^diag{-l,+l,+l). 

The Christoffel symbol corresponding to this flat "supermetric" vanishes, V^q — 0, and 
the Hamiltonian equations translate into 

d\^ 2 ^dg^' ' 

I.e. it is a kind of Newton'second law with non-vanishing force term appearing on the 
right hand side, d'^n/dX^ = -l&JZ/dQ, d'^p+fdX'^ = |97^/9/3+, d^p^/dX^ = |97^/9/3_. 

However, the configuration space variables g^ = (n,/9+,/3_) may be mapped to a 
(affinely parametrized) geodesic flow on the configuration space manifold equipped with 
the (mini)superspace metric 

Gab = TZGab = -e^"(V^(/3+, /?_) - l)diag{-l, 1, 1) (3.17) 

with the parametrization 

dX = ndX = -e^^{V{p+, /?_) - l)dX (3.18) 
along the trajectory of the mixmaster collapse. Thus, we have 

^ + f^^^^ = (3.19) 
dX dX ^ ^ 
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This procedure introduces a host of singularities at points where 7Z = i.e. where V = 1 
(see also an interesting discussion in Misner pp. 453-454). 



The fact, recently emphasized also by Burd and Tavakol that the singularities 
introduced by the supermetric Gijki prevents us from mapping the entire mixmaster col- 
lapse orbit into one single, unbroken and simple geodesic orbit - is of course an important 
obstacle to this approach and was, in fact, previously emphasized by one of us (S.E.R.). 
Despite these troubles, let us suppose that we nevertheless carry out these transformations 
at points in the mixmaster configuration space = (fi, where the superspace 
metric Gab = Gijki do not introduce artificial singularities, i.e. at those points where 
V{P+,P-) 7^ 1. With this supermetric Gijki one may thus construct non-trivial intervals 
in which the mixmaster metric is mapped into geodesies (in Wheeler's (mini) superspace). 
We shall start with some remarks on the properties of the distance measure between two 
nearby mixmaster metrics ^^^gix and ^^^gix + S^^^gjx- Let us first note, that the distance 
between g"^ = {Q, and g"^ + 6g^ = {Q + 6Q, (3+ + 5(3+, (3- + 6(3-) is of Lorentzian 

signature. This indefiniteness is a property of the supermetric, both in its original form 

I ISg'^l |2 = ds^ = GABSg^Sg^ = -dQ^ + d(3l + d(3l (3.20) 

and in the conformally rescaled form 

||5^^||2 = d~s^ = GABSg^Sg^ = -e~^^{V{(3+, (3^) - l){-dQ^ + d(3l + d(3l). (3.21) 

Thus the distance can take positive, zero or negative values and one has the unfamiliar 
situation (relative to positive definite Euclidean distance measures implemented in the 
context of non-relativistic Hamiltonian dynamical systems) - as previously stressed - that 
two configuration space points (mixmaster three- metrics) which is not identical may have 
zero distance with respect to these naturally induced distance measures. 

These distance measures have on the other hand the good property of being invariant 
under canonical coordinate transformations. Thus, if we have a change of coordinates 
(configuration space variables), g^ — > g*^ the distance measure ds^ = GAB^g^^^g^ is 
invariant 

\\6g^\\' = ds' = GAB^gHg^ = G\j,5g*Hg*^ = \\6g*^\\' 
under such transformations, since Gab transforms properly as a tensor. 

One can easily see this by recalling that the coordinate transformation g"^ g*^ 
induces the canonical transformation of momenta 



Pa^Pa = tttj Pb 



dg^ 
dg* 



Therefore the Hamiltonian Ti reads 

n = -^{G^^'PAPB -n) = -^{G^^'i^-I^PcPD -n) = -^{G^^^'p^fs - 

thus justifying claims of the formula ( p.22|) . 
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A quantity which is invariant under a large class of coordinate reparametrizations 
(canonical transformations in the Hamiltonian formulation) was considered in works by 



Biesiada, Lapeta, Szcz§sny and Szydlowski 



where the Ricci scalar of the 
manifold on which the mixmaster model acts as a geodesic flow has been extracted. As 
we have seen here this is precisely to extract the Ricci scalar for the Wheeler's superspace 
metric Gab (in the conformally rescaled version where it describes the three-metrics as 
geodesies) and negative values R < are naturally interpreted as a local measure of 
exponential instability via the geodesic deviation equation with respect to the conformally 
rescaled distance measure - of Lorentzian signature - on the space of mixmaster three- 
metrics. 

Despite the problems with unwanted singularities induced by this approach (as noted 
by e.g. Misner |2^) we believe that insufficient attention has been paid towards the 
use of distance measures naturally induced by the structure of general relativity, e.g. the 
Wheeler's superspace metric as a distance measure on the space of mixmaster collapses (in 
the context of, for example, discussing the chaotic properties) - instead of using artificial 
and completely arbitrary Euclidean distance measures |P, |10], Q on the solution space 
(which are not supported by the structure of general relativity). 

As concerns our discussion of the mixmaster gravitational collapse, the (mini) superspace 
Ricci scalar R may in principle be calculated directly from the conformally rescaled 
Wheeler's (mini)superspace metric, 
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However, ADM variables are not very suitable for such calculations due to the complicated 
form of the potential term. Since our aim is to extract the Ricci scalar R{Maupertuis), 
which is an invariant under canonical coordinate transformations, it should not matter in 
which representation we work. The ADM variables and the Bogoyavlenskii variables are 
related by canonical coordinate transformations.^ 

Therefore, as in pO|, we rather consider the Bogoyavlenskii Hamiltonian in (a, b, c) 
variables 
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"^{PaPha^b'^ + PbPcb'^C^ + PaPcO^C^^ 
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(for the mixmaster collapse) in which case the potential is less complicated at the expense 
that the metric read off from the kinetic energy is no longer diagonal, 
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More precisely, one may transform the ADM variables to a,/3,7 variables (cf. equation (3.14)) and 
make a corresponding canonical coordinate transformation of the ADM momentum variables to , p/3 , . 
These are then the Bogoyavlenskii momentum variables, cf. Bogoyavlenskii [E6|, p. 40, and one regains 



the Bogoyavlenskii Hamiltonian from the ADM Hamiltonian up to a conformal factor a^b^c^ 
which is not important because the Hamiltonians are equal to zero, H — 0. 



„2(a+/3+7) 
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Explicit evaluation of the Ricci scalar ( p.l2| ), corresponding to the conformally rescaled 
metric, induced by the kinetic term in the Bogoyavlenski Hamiltonian, yields 



Rj^ = -^—^{a^ + 6« + c« - 2a%^ - 2a^c^ - 2h^c^ + IQa^^^a^ + 6^ + ^)) (3.23) 

8(— K j'^ 

where 

-V = +\g{^R) = -^{a^ + h^ + c^- 2a%^ - 2a''c^ - 21? c^) (3.24) 

As concerns the aforementioned instability criterion applied to the mixmaster three- 
geometries the idea is now to find regions in the configuration space where i? < and see 
whether a typical trajectory is "confined" to such regions or just traverse them quickly. 
Clearly, R cannot be negative all over the entire minisuperspace. For example, it is easy 
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to see that i? > in the case of isotropy, a = b = c. On the other hand one observes pO 
that -R < in the asymptotical BKL regime, i.e. when a » b,c (or cyclic permutations 
thereof). 

For sake of illustration we performed a numerical experiment in which we have inte- 
grated the vacuum mixmaster field equations 



(3.25) 



supplemented by a first integral constraint 

/ = arPr + drlr + Mr - ^(fl^ + &^ + C^) - 2{a^h'^ - 2o? - 2\? C^) = . 

Here t = J dt/abc denotes the standard time variable, and subscript r means differ- 
entiation with respect to the time variable r. We used a numerical code from S.E. Rugh 
[0 based on the fourth-order Runge-Kutta integrator with check for the first integral 
constraint (for details see ||12|| ). 

The figure displays the temporal (in r-time) behavior of the minisuperspace Ricci 
scalar as felt by a trial trajectory. We have selected a set of reference initial conditions 
as in A. Zardecki but have adjusted the value of c' to make the first integral vanish 
to machine precision. (Such an adjustment is indeed necessary. Cf. discussions in S.E. 



Rugh [|T2| and D. Hobill [|13])- This yields the starting conditions 
a = 1.85400.. b = 0.438500.. , c = 0.085400.. 



a 



' - 0.429200.. , b' = 0.135500.. , c = 2.964843279 (3.26) 



On the figure the minisuperspace Ricci scalar is displayed. For comparison also the 
evolution of (logarithmic) scale factors a, (3, 7 is superimposed. It turns out that R is 
negative all the time but varies over many orders of magnitude. The latter is connected 
with the singularities generated by a vanishing conformal factor at points where the 
potential is zero. At these points R goes to (minus) infinity and a numerical cure against 
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it was to provide an arbitrary cut off. This means that our trial trajectory of the mixmaster 
collapse stays in the regions of minisuperspace where "local instability" occurs. Of course 
such an illustration cannot be representative for any larger set of trajectories. For example, 
the minisuperspace Ricci scalar calculated along a Taub axisymmetric solution was found 
to be positive. However, even in cases when i? < it is not possible to claim, cf. our 
previous discussion about local and global properties, that this fact accounts for global 
chaotic behavior - even though some may be tempted to make such an inference. 



4 Concluding remarks. 

The possibility of mapping the dynamics of a wide class of Hamiltonian systems to geodesic 
flows on Riemannian manifolds by application of Maupertuis principle, and extracting 
coordinate invariant information about local instability of the trajectories (from the cor- 
responding Riemann tensor) has been applied in various contexts, ranging from N-body 
systems (a gas) interacting via Van der Waals forces (Krylov |2^), the instability (relax- 
ation) properties of a collisionless gas interacting via gravitational forces (Gurzadyan and 
Savvidy [^) and N-body systems interacting via the Debye-Hiickel potential given by 
V{r) ~ e^'^'^/r (a model for a hot dilute plasma. Van Velsen [p5|] ). 

In this paper we have examined virtues and drawbacks of applying this idea to general 
relativity |18[-||22[. Coordinate invariant information about local instability properties is 



in this case - as in the applications to the aforementioned non-relativistic Hamiltonian 
systems - imagined to be extracted either from negativity of the Ricci scalar (which is 
the sum over principal sectional curvatures) or by considering the (principal) sectional 
curvatures themselves |^ . 



We have shown that the Maupertuis principle as a way of geometrizing the Hamiltonian 
dynamics when implemented to general relativity is actually reproducing Wheeler's idea of 
a superspace in which three-geometries evolve along geodesic lines. We have examined this 
superspace in order to get an insight into the construction of gauge invariant (geometrical) 
criteria for local instability in general relativity. The original dream of having Einstein's 
equations acting as a geodesic flow has been abandoned by its inventors a long time 
ago because of artificial singularities created by the conformal rescaling of the metric in 
points where the ADM potential equals zero. In so far as for its original purpose i.e. 
to provide a suitable state-space for quantization of gravity, people could (and actually 
did) retreat to more secure ground of the weaker form of the superspace (giving up the 
geodesic picture) the conformally resettled superspace is crucial for our proposal. The idea 
is simple: If one can cast the dynamics into a geodesic flow then the sensitive dependence 
on initial conditions (a key ingredient of instability and chaos) may be captured at the 
level of Jacobi equation for geodesic deviation. The original intention was to point to 
a certain reasonable approach rather than provide an ultimate solution (if any exist) to 
the problem of invariant description of chaos. We have noticed some obvious dangers 
when one makes global conclusions from local considerations. This issue is far from being 
trivial however, since the history of differential geometry teaches that starting from Gauss' 
theoremtt egregium a lot of effort (sometimes spectacularly successful) has been paid to 
the question of how to derive global characteristics from local ones. 

We have shown that expressions - used previously - such like (|2.22|) relating directly the 
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Ricci scalar with Lyapunov exponents ||T9|, ^ are not true. Therefore, unfortunately, 
one cannot base a short time average (STA) Lyapunov exponent on this expression either, 
cf. Burd and Tavakol The STA approach as a tool for studying non-uniform 

dynamics is based on evaluation of the Lyapunov exponents (extracted in a standard way 
from the Jacobian) over a finite time interval. The mixmaster toy-model gravitational col- 
lapse examined here as an example of a spacetime metric with complex (chaotic) behavior 
is in fact rather well understood — the BKL-description is derived under assumptions 
which are very good and indeed numerically confirmed in Rugh [|T2| and by Berger . 
Nevertheless, if we want to characterize - or rather to develop measures of - chaos in 
the highly non-linear Einstein equations, the mixmaster collapse is a suitable toy-model 
laboratory to address such questions: If we are not even able to invent useful indicators 
of chaos in the moderately simple example of the mixmaster space-time metric, how are 
we ever going to be able to deal with these issues for more general and more complicated 
space-time metrics? 

The observation of gauge variance of standard indicators like that of a Lyapunov 



exponent in the context of general relativity was emphasized by Rugh |12|, |1J] and also 
by Pullin |jl5[. (See also recent discussion in Rugh fl^). If we apply standard techniques 
of extracting the Lyapunov exponents we at first sight have chaos in some gauges and 
in other gauges there seems to be no chaos. Apparently, this situation arises because 
we are looking at the problem in the wrong way - we have not posed a "gauge invariant 
question" ! 

Whereas Rugh |14] was expressing hope of constructing indicators which are invariant 
under a large class of coordinate transformations, i.e. indicators which do not refer to a 
"particular gauge" (a problem which is easier to point out than to solve) it appears that 
Pullin |]T3[ was resorting to the Poincare disc as a kind of "selected set of coordinates" 
— a similar viewpoint is shared by Misner ||4^. It appears to us that no "gauge" is 
better than others (is chaos in general relativity a concept which should only be defined 
in certain selected frames of coordinates?) - yet some gauges may make things partic- 
ularly simple or elegant thereby expressing the explicit or hidden symmetry (geometric 
structure) of the problem. The geometric approach (via Maupertuis principle) discussed 
in the present paper also points toward a certain structure uncovered from the Hamil- 
tonian formulation of general relativity (which turned out to be equivalent to Wheeler's 
superspace) as a "natural" one. The problem is however not with "preferred" gauges but 
rather whether the criteria used to establish certain properties (such like chaoticity) are 
gauge-invariant themselves. Therefore we ought to seek measures of chaos (if they are 
possible to construct) which are invariant under some large class of gauge transformations. 

Here we have examined an approach which points to a more gauge invariant con- 
struction involving Wheeler's superpace. However, this approach does (also) have several 
obstacles, e.g. when implemented in the toy-model context of the mixmaster gravitational 
collapse. Hence it is safe to conclude that how to characterize chaos most elegantly in 
general relativity (even in the context of simple toy-models) remains an open question. 
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Figure 1: An illustration of the mixmaster minisuper space Ricci scalar for a trial trajectory. 

Flat minima denote that the depth is out of range of the figure. The Ricci scalar is negative 
all the time (but rapid changes over many orders of magnitude makes it difficult to perceive 
this fact on the figure). Indeed the artificial singularities introduced by the conformally rescaled 
superspace metric (see the text) make the Ricci scalar i? go to — oo while traversing the null- 
surface on which the potential is zero. 



Figure caption 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9408030vl 



